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'sj" ■ Abstract. Ideal hydroynamics provides an excellent description of all aspects of the single-particle 

I spectra of all hadrons with transverse momenta below about 1 .5-2 Ge V/c at RHIC. This is shown to 

i require rapid local thermalization at a time scale below I fm/c and at energy densities which exceed 

the critical value for color deconfinement by an order of magnitude. The only known thermalized 
state at such energy densities is the quark-gluon plasma (QGP). The rapid thermalization indicates 
that the QGP is a strongly interacting liquid rather than the weakly interacting gas of quarks and 
gluons that was previously expected. 

(N 

<N ■ Collective flow as an indicator of thermalization 

> ' 

\0 • Transverse flow in heavy-ion collisions is an unavoidable consequence of thermal- 

O I ization. Thermalization generates thermodynamic pressure in the matter created in the 

' collision, which acts against the surrounding vacuum and causes rapid collective expan- 

. sion ("flow") of the reaction zone. Since the quark-gluon plasma (QGP) is a thermalized 

system of deconfined quarks, antiquarks, and gluons, collective flow is a necessary result 
of QGP formation in heavy-ion collisions, and its absence could be taken as proof that 
no such plasma was ever formed. Its presence, however, does not automatically signal 
O I QGP formation. Detailed studies of the observed final state flow pattern are necessary 

^ I to convince oneself that the reflected time-integrated pressure history of the collision re- 

gion indeed requires a thermalized state in the early collision stage whose pressure and 
energy density are so high that it can no longer be mistaken as consisting of conventional 
^ ' hadronic matter. 

c3 ■ Whereas radial flow (the azimuthally symmetric component of collective expansion 

transverse to the beam direction) integrates over the entire pressure history of the ex- 
panding fireball, anisotropic elliptic flow ^ (quantified by the second harmonic co- 
efficient V2{yjP±',b) of a Fourier expansion in of the measured hadron spectrum 
dN/{dyp \ dp±d^p) [2]) is strongly weighted towards the very early stages of the ex- 
pansion [3]. The higher the initial energy density, the less contributions it receives from 
the late, hadronic stage of the collision. At RHIC energies, elliptic flow almost satu- 
rates |i4J before the energy density has dropped to the critical value Ccr ~ 1 GeV/fm-^ 
(Tcr ~ 170 MeV) where normal hadrons can begin to form |5]. At higher LHC ener- 
gies, the elliptic flow is expected to peak even before the onset of hadronization (see 
curve c in Fig. 7 of f?]). At sufficiently high collision energies, elliptic flow is thus 
a QGP signature, probing the QGP equation of state p{e). Its advantage over other 
"early signatures" is that it affects the bulk of the hadrons and can thus be measured 
differentially with high statistical accuracy. 



The reason why elliptic flow must develop early in the collision is easy to understand. 
Since individual nucleon-nucleon collisions produce azimuthally symmetric momentum 
spectra, any final state momentum anisotropics must be generated dynamically during 
the nuclear reaction. They require the existence of an initial spatial anisotropy of the 
reaction zone, either by colliding deformed nuclei such as U+U or by colliding 

spherical nuclei at non-zero impact parameter b^^O (the practical method of choice 
so far). Final state interactions within the produced matter transfer the initial spatial 
anisotropy onto a final momentum anisotropy. Microscopic transport calculations 
show a monotonic dependence of V2 on the opacity (density times scattering cross 
section) of the produced matter which is inversely related to its thermalization time. 
These studies strongly suggest that, for a given initial spatial anisotropy fix, the maximum 
momentum-space response V2 is obtained in the ideal hydrodynamic limit which assumes 
perfect local thermal equilibrium at every space-time point (i.e. a thermalization time 
which is much shorter than any macroscopic time scale in the system). Any significant 
delay of thermalization (modelled, for example, as an initial free- streaming stage) causes 
a decrease of the initial spatial anisotropy without concurrent build-up of momentum 
anisotropics, thereby reducing the finally observed elliptic flow signal [4]. 

In this talk I will present results from hydrodynamic simulations of hadronic spectra 
and elliptic flow at RHIC energies. We will see that the hydrodynamic approach pro- 
vides an excellent quantitative description of the bulk of the data and fails only for 
very peripheral Au-i-Au collisions and/or at high pj^> 1.5—2 GeV/c. That the hydro- 
dynamic approach fails if the initial nuclear overlap region becomes too small or the 
transverse momentum of the measured hadrons becomes too large is not unexpected. 
However, where exactly hydrodynamics begins to break down gives important infor- 
mation about the microscopic rescattering dynamics. What is really surprising is that 
the hydrodynamic approach works so well in central and semi-central collisions where 
it quantitatively reproduces the momenta of more than 99% of the particles. Below 
p^ = l.5 GeV/c the ellintic flow data lEi E Q actually exhaust the hydrodynami- 
cally predicted [l^ [l^ [la, upper limit. The significance of this agreement can 
hardly be overstressed, since it implies one of the biggest surprises so far at RHIC: the 
produced matter (which I call QGP since this is only known viable concept of thermal- 
ized matter at e ~ (10—20) ecr) is not the originally expected weakly interacting gas of 
quarks and gluons (wQGP for "weakly interacting quark-gluon plasma"), but a strongly 
coupled liquid with extremely small viscosity (sQGP for "strongly interacting quark- 
gluon plasma" [18]). In fact, upper limits [19] on the dimensionless ratio of viscosity to 
entropy density, rj/s, based on an analysis of the RHIC elliptic flow data, indicate that 
the sQGP is less viscous, by about an order of magnitude, than even liquid helium below 
the transition to superfluidity. The sQGP is therefore the most ideal fluid ever observedl 



Hydrodynamic expansion in heavy-ion collisions 

The natural language for describing collective flow phenomena is hydrodynamics. Its 
equations control the space-time evolution of the pressure, energy and particle densities 
and of the local fluid velocity. The system of hydrodynamic equations is closed by 



specifying an equation of state which gives the pressure as a function of the energy and 
particle densities. In the ideal fluid (non-viscous) limit, the approach assumes that the 
microscopic momentum distribution is thermal at every point in space and time (note that 
this does not require chemical equilibrium - chemically non-equilibrated situations can 
be treated by introducing into the equation of state non-equilibrium chemical potentials 
for each particle species L20i. .2L .2 2 ]). Small deviations from local thermal equilibrium 
can in principle be dealt with by including viscosity, heat conduction and diffusion 
effects, but such a program is made difficult in practice by a number of technical and 
conceptual questions [12311 and has so far not been successfully applied to relativistic 
fluids. Stronger deviations from local thermal equilibrium require a microscopic phase- 
space approach (kinetic transport theory), but in this case the concepts of equation 
of state and local fluid velocity field themselves become ambiguous, and the direct 
connection between flow observables and the equation of state of the expanding matter 
is lost. 

The assumption of local thermal equilibrium in hydrodynamics is an external input, 
and hydrodynamics offers no direct insights about the equilibration mechanisms. It is 
clearly invalid during the initial particle production and early recattering stage, and it 
again breaks down towards the end when the matter has become so dilute that rescat- 
tering ceases and the hadrons "freeze out". The hydrodynamic approach thus requires 
a set of initial conditions for the dynamic variables at the earliest time at which the as- 
sumption of local thermal equilibrium is applicable, and a "freeze-out prescription " at 
the end. For the latter we use the Cooper-Frye algorithm [24] which implements an ide- 
alized sudden transition from perfect local thermal equilibrium to free- streaming. This 
is not unreasonable because freeze-out (of particle species /) is controlled by a compe- 
tition between the local expansion rate d ■u{x) (where u^{x) is the fluid velocity field) 
and the local scattering rate Y.j{OijVij)pj{x) (where the sum goes over all particle species 
with densities pj{x) and {OijVij) is the momentum-averaged transport cross section for 
scattering between particle species i and j, weighted with their relative velocity); while 
the local expansion rate turns out to have a rather weak time-dependence, the scattering 
rate drops very steeply as a function of time, due to the rapid dilution of the parti- 
cle densities pj |25J, causing a rapid transition to free-streaming. - A better algorithm 
lE^ I14I1 switches from a hydrodynamic description to a microscopic hadron cascade at 
or shortly after the quark-hadron transition, before the matter becomes too dilute, and 
lets the cascade handle the freeze-out kinetics. This also correctly reproduces the final 
chemical composition of the fireball, since the particle abundances already freeze out at 
hadronization, due to a lack of particle-number changing inelastic rescattering processes 
in the hadronic phase. The resulting radial flow patterns 1I14I] from such an improved 
freeze-out algorithm don't differ much from our simpler Cooper-Frye based approach. 



Hydrodynamic radial flow and RHIC particle spectra 

We have solved the relativistic equations for ideal hydrodynamics, as described in 
10]. To simplify the numerical task, we imposed boost-invariant longitudinal expansion 
analytically |[il0]. As long as we focus on the transverse expansion dynamics near 



midrapidity (the region which most RHIC experiments cover best), this does not give up 
any essential physics. 

We use an EOS which is constructed by matching a free ideal quark-gluon gas above 
Tcr to a realistic hadron resonance gas [28, 29] below T^r, using a Maxwell construc- 
tion and fixing Tcr = 165 MeV to reproduce lattice QCD results [5]. The Maxwell con- 
struction leads to an artificial first order transition, with a latent heat of 1.15 GeV/fm^, 
whereas the lattice QCD data indicate either a very weak first order transition or a rapid, 
but smooth, crossover. As long as this crossover is very steep (as the lattice data indi- 
cate Jil]), the dynamical consequences of our first-order idealization of the EOS are not 
expected to be significant. 

We have used two variants of the hadron resonance gas below Tex'. In our first sets of 
calculations 1999-2002, we assumed the hadron resonance gas to be in chemical equi- 
librium until kinetic freeze-out. RHIC data on particle abundance ratios show, however, 
that hadronic particle yields freeze out directly at T^x 1 30], due to inefficiency of particle- 
changing inelastic reactions in the relatively dilute and rapidly expanding hadron gas 
phase below T^r. More recent hydrodynamic calculations ll2lL l3lll therefore include non- 
equilibrium chemical potentials ll20ll2ll '221 ensuring number conservation for individ- 
ual (stable) hadron species below T^x- This does not affect the equation of state pie ) 
and therefore leads to the same flow pattern as the chemical equilibrium EOS 11211]. 
What does change, however, is the relationship between energy density e and temper- 
ature T since in the chemical non-equilibrium case more of the energy is stored in the 
rest masses of heavy baryons and antibaryons (which are not allowed to annihilate as the 
temperature drops). The same decoupling energy density e^ec = 0.075 GeV/fm^ obtained 
from a fit to RHIC data (see next paragraph) therefore corresponds to r^ec = 130 MeV 
in the chemical equilibrium case, but drops to the much lower value Tjec = 100 MeV if 
hadronic chemical freeze-out at = 165 MeV is properly taken into account li2Li3l1l . 

The initial and final conditions for the hydrodynamic evolution are fixed by fitting 
the total charged multiplicity dNch/dy and the pion and proton spectra at midrapidity 
in central {b = 0) collisions. (All calculated hadron spectra include feeddown from de- 
cays of unstable hadron resonances [32].) Since these two hadrons have quite different 
masses, their spectra allow for the independent extraction of the temperature and aver- 
age radial flow velocity at freeze-out 113 3ll . The freeze-out temperature together with the 
total charged multiplicity determines the total fireball entropy which in ideal hydrody- 
namics is preserved during the expansion. During the early, predominantly longitudinal 
expansion stage, where the fireball volume increases linearly with proper time T, this 
constraint fixes the product T ■ / d^r±s{r±, t) where s{r±, t) is the entropy density dis- 
tribution in the plane transverse to the beam. Since the initial shape of this distribution is 
fixed by the nuclear overlap geometry, using a Glauber model ansatz [ 17], this constraint 
determines the normalization of s{r±, Teq) (i.e. the central entropy density 50 = ^(0, Teq)) 
as a function of the initial thermalization time Teq- The Glauber model contains one 
additional parameter: the ratio of "soft" and "hard" components of initial particle pro- 
duction. Since these scale with the transverse density of wounded nucleons and binary 
nucleon-nucleon collisions, respectively, their ratio can be determined from the collision 
centrality dependence of the produced charged multiplicity dN^h / dv aim . We use 25% 
hard and 75% soft contributions to the initial entropy production ll34ll . The initial ther- 
malization time Teq, finally, is fixed by the need for the hydrodynamic evolution to have 



enough time to generate the finally observed radial flow. At RHIC energies, we find that 
this requires thermalization times Teq<0.6fm/c. In our calculations we take the upper 
limit of this interval. In Au+Au collisions at v^= 100 A GeV, the initial central entropy 
density at this time is 5o= 1 10 fm^-^, corresponding to an initial central energy density' 
eo = 30 GeV/fm^ ^ 30 and an initial central temperature Tq = 360 MeV ^ITcr-^ 
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FIGURE 1. (color online) Negative pion, kaon, antiproton, and Q. spectra from central Au+Au colli- 
sions at ^/i = 200 A GeV, as measured by the four RHIC experiments I35ll36ll37l 138113911 . The curves show 
hydrodynamical calculations as described in the text. 

In Fig.[T]l show the (absolutely normalized) single particle p^-spectra for negatively 
charged pions, kaons and antiprotons (left panel) as well as Q. baryons (right panel) 
measured in Au+Au collisions at RHIC together with hydrodynamical results |31]. In 
order to illustrate the effect of additional radial flow generated by elastic scattering in 
the late hadronic stage below T^, two sets of curves are shown: the lower (blue) bands 
correspond to kinetic decoupling at 165 MeV, whereas the upper (red) bands as- 
sume decoupling at Tdec = 100 MeV. The width of the bands indicates the sensitivity of 
the calculated spectra to an initial transverse flow of the fireball already at the time of 
thermalization: The lower end assumes no initial transverse flow whereas the upper end 
implements an initial radial flow profile Vr{r±, Teq) = tanh(arj^) with a = 0.02fm^^ at 
T'eq = 0.6fm/c (which seems to be slightly preferred by the data). The hydrodynamic 
model output shows ||4j] that it takes about 9-10 fm/c until the fireball has become suf- 
ficiently dilute to completely convert to hadronic matter and another 7-8 fm/c to com- 
pletely decouple. Figured] shows clearly that by the time of hadronization hydrodynam- 
ics has not yet generated enough radial flow to reproduce the measured proton and Q. 
spectra; these heavy hadrons, which are particularly sensitive to radial flow effects, re- 
quire the additional collective "push" created by resonant (quasi)elastic interactions dur- 
ing the fairly long-lived hadronic rescattering stage between 7^ and Tdec- Even though 
independent flow fits to Tt,K,p spectra on the one hand and multistrange hyperon spec- 
tra on the other [40] suggest (within some systematic uncertainty related to the strong 



' When averaged over the transverse profile, this corresponds to (e) = 13GeV/fm^ which is still a good 
order of magnitude above the critical value for deconfinement, ga- = 0.6 — 1 GeV/fm^ j^l • 
^ If the initial QGP is not chemically equilibrated, but rather dominated by gluons, this initial temperature 
could be as high as 460 MeV. 



anticorrelation between the extracted flow and temperature values) that S and Q. hyper- 
ons decouple slightly earlier (at somewhat higher temperature and with less radial flow) 
than pions, kaons, and protons, their immediate decoupling directly at hadronization is 
obviously not dynamically consistent with the successful hydrodynamic approach. 

The strong flattening of the (anti)proton spectra by radial flow provides a natural ex- 
planation for the (initially puzzling) experimental observation that for p±>2 GeV/c an- 
tiprotons become more abundant than pions. For a hydrodynamically expanding ther- 
malized fireball, at relativistic transverse momenta p± ^ mo all hadron spectra have 
the same slope ll33ll . and at fixed m^^mo their relative normalization is given by 
{giXi) / {gjXj) (where gi_j is the spin-isospin degeneracy factor and A/.j = e^'j/^ is the 
fugacity of hadron species ?, j). Due to the large antiproton chemical potential at kinetic 
freeze-out which is required to maintain the p abundance at temperatures below chem- 
ical freeze-out, this asymptotic p/iz^ ratio is predicted to be about 28 |41]. Within the 
hydrodynamic approach the surprising fact is therefore not that p/n>la.ip>^^>2 GeV/c, 
but that this ratio seems to saturate around 1 and never much exceeds this value llssll . 
As we will see below, this is a signature of the beginning breakdown of the hydrody- 
namic model which stops working for baryons with transverse momenta above about 
2.5 GeV/c. 

As shown elsewhere (see Fig. 1 in lH^]), once the hydrodynamic model parameters 
have been fixed to describe pions, protons, and total multiplicity in central Au-i-Au 
collisions, the model describes these and all other hadron spectra not only in central, 
but also in peripheral collisions, up to impact parameters of about 10 fm and with 
similar quality. No additional parameters enter at non-zero impact parameter - only the 
initial conditions change due to the changing overlap geometry, but this is completely 
accounted for by the Glauber model. Therefore, hydrodynamic results for the elliptic 
flow, discussed in the next Section, are parameter-free predictions of the model. 

Hydrodynamic elliptic flow and RHIC data 

Figure El shows the predictions for the elliptic flow coefficient V2 from Au-i-Au col- 
lisions at RHIC, together with the data [11, ij,!!^. For impact parameters b<lim 
(corresponding to Rch/ '^max > 0-5) and transverse momenta p^<l.5 — 2 GeV/c the data 
are seen to exhaust the upper limit for V2 obtained from the hydrodynamic calculations. 
For larger impact parameters b>lim the p^-averaged elliptic flow V2 increasingly lags 
behind the hydrodynamic prediction. It is tempting to attribute this to a lack of early 
thermalization when the initial overlap region becomes too small ll34ll . However, the 
work by Teaney Hill , who coupled hydrodynamic evolution of the quark-gluon plasma 
above T^t with a microscopic kinetic evolution of a hadronic resonance gas using RQMD 
below Tcr, showed that at least a strong contributing factor to the lack of elliptic flow in 
peripheral collisions (as well as in central collisions at lower energy) is the large viscos- 
ity in the late hadronic stage. If the initial energy density is not high enough, the initial 
spatial deformation does not completely disappear until hadronization. Whereas ideal 
hydrodynamics then predicts a continued growth of the elliptic flow during the hadronic 
stage, driven by the still existing anisotropics in the pressure gradients, realistic hadron 
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FIGURE 2. (color online) Elliptic flow coefficient V2 for all charged particles from 130 A GeV Au+Au 
collisions (left panel (v^) and for several identified hadron species from 200AGeV Au+Au collisions 
(right panel |44]), compared with hydrodynamic predictions. The left panel shows the /?x-averaged 
elliptic flow as a function of collision centrality, parametrized by the charged multiplicity density «ch 
at midrapidity («max corresponds to the largest value in central collisions). The right panel shows the 
differential elliptic flow V2{p±) for minimum bias collisions. The data were collected by STAR and 
PHENIX llJlHlil. The curves in the left panel are hydrodynamic calculations corresponding to 
different choices for the initial energy density profile (see 1 17] for details). The curves in the right panel 
were first published in fl^ . 



cascades produce very little additional elliptic flow, indicating that the hadronic matter is 
very viscous. To completely saturate the elliptic flow value predicted by ideal hydrody- 
namics thus requires not only very rapid thermalization at the beginning of the reaction, 
but also a sufficiently long lifetime of the thermalized QGP stage such that the initial spa- 
tial anisotropy has time to fully decay and the momentum anisotropy can fully develop 
before the system enters the viscous hadron gas phase. This explains why the elliptic 
flow lags behind the hydrodynamic prediction at RHIC in peripheral Au-i-Au collisions 
at midrapidity, in minimum bias Au-i-Au collisions at forward and backward rapidities 
lEllE5l t46]. as well as quite generally at lower collision energies |S3.Ss|]- From an anal- 
ysis of all the available data we recently concluded [46] that an average initial energy 
density {e)x^^ > 10 GeV/fm^ appears to be necessary to provide sufficient evolution time 
in the plasma phase for full development of the momentum anisotropy before hadroniza- 
tion and therefore full saturation of the hydrodynamic elliptic flow limit. 

The -dependence of V2 (right panel of Fig. |2l) tells another story: whereas the 
hydrodynamic model predicts a continuous rise of V2 with increasing p^, the measured 
elliptic flow appears to saturate at high p±. Baryon elliptic flow (here shown are protons 
and A hyperons, but E and Q. hyperons follow the same systematics [l49tl ) saturates at a 
higher value than meson (pion and kaon) elliptic flow. These saturation patterns imply a 
break of the data away from the hydrodynamic prediction, and again this break happens 
at higher for baryons {p^l^^ ~ 2.3 - 2.5 GeV/c) than for mesons (p^^^^^ ~ 1 .5 GeV/c). 
Below these break points, the elliptic flow of all hadrons measured so far is very 
well described by hydrodynamics. In particular the hydrodynamically predicted mass 
splitting of V2 at low p± is perfectly reproduced by the data. Hydrodynamics thus 
gives an excellent description of all hadron spectra below p± = l.5 GeV/c. Due to the 



exponentially falling spectra, this includes more than 99% of all produced hadrons, 
so that it is fair to say that the bulk of the fireball matter formed in Au+Au collisions at 
RHIC behaves like a perfect fiuid. 

In jh we showed that the slope of V2{pi_) and its mass splitting are sensitive to 
the equation of state. The mass splitting is larger for an EOS which includes a quark- 
hadron phase transition than for a pure resonance gas EOS extrapolated to arbitrarily 
high temperatures. Although requiring many more systematic model studies than so 
far available, these features offer the perspective of extracting detailed knowledge on 
the EOS during the early expansion stages from combined radial and elliptic flow 
measurements. So far we only know qualitatively (see Fig. 2 in [42] and the discussion 
in fl4]) that, for given angle-averaged p^ spectra, the measured elliptic flow prefers an 
EOS that includes a phase transition at T^x = 165 MeV over one without such a transition. 

The excellent agreement between data and hydrodynamic model at transverse mo- 
menta below about 2 GeV becomes even more impressive after one begins to realize how 
easily it is destroyed: Parton cascade simulations with standard HIJING input generate 
almost no elliptic flow and require an artificial increase of the opacity of the partonic 
matter by a factor 80 to reproduce the RHIC data [9]. This shows that the QGP is much 
more strongly coupled than achievable within a standard perturbative QCD approach. 
Hadronic cascades of the RQMD and URQMD type (in which the high-density initial 
state is parametrized by non-interacting, pressureless QCD strings) predict Iscill too little 
elliptic flow and a decrease of V2 from SPS to RHIC, contrary to the data. Generically, 
models that don't include very strong rescattering in the very early fireball expansion 
stage fail to reproduce the elliptic flow data. 

As I will show in the next section, any significant delay of thermalization immediately 
leads to a significant loss of elliptic flow signal - a fact that can be used to establish 
tight upper limits on the thermalization time scale. Even partial thermalization doesn't 
work: In Isill we studied the transverse hydrodynamic evolution of a toy model in 
which only transverse momenta were thermalized while in longitudinal direction the 
fireball evolved ballistically. (One can think of this as implementing into the transverse 
hydrodynamics a very large value for the shear viscosity.) Since in this model there is 
no longitudinal pressure, more work is done by the pressure in the transverse direction, 
and one has a stiff er effective equation of state p = e/2 (instead of the ideal quark-gluon 
gas EOS p — e/3) driving the transverse expansion. To obtain the same pj^-spectra as 
before one must thus retune the initial conditions and let the hydrodynamic expansion 
start later, in order to avoid getting too much radial flow. The consequence of this is a 
reduction of the elliptic flow V2 by almost a factor 2, to levels far below the experimental 
data. This exercise demonstrates a strong sensitivity of the elliptic flow to early shear 
viscosity. A complementary analysis by Teaney, who discussed late-stage viscous effects 
on hadron spectra, elliptic flow and HBT radii in [19], arrives at similar conclusions. This 
sensitivity of V2 to shear viscosity opens the perspective of limiting the shear viscosity 
of a quark-gluon plasma experimentally by performing systematic comparisons between 
elliptic flow data and viscous (i.e. non-ideal) hydrodynamics. This has recently led to 
renewed interest in practical implementations of viscous relativistic fluid dynamics ll52ll . 



Upper limits for the QGP thermalization time 



The agreement between the elliptic flow data at RHIC and the hydrodynamic simula- 
tions can be used to set an upper limit on the thermalization time scale in the dense par- 
tonic matter formed in the collision which is independent from the value Teq < 0.6 fm/c 
obtained two chapters earlier by fitting the model to central collision spectra. Let me 
illustrate this by referring to Fig. |3] which, to some extent, is a different way of plotting 
the results shown in the left panel of Fig. |2l In ideal hydrodynamics, the final elliptic 
flow V2 is directly proportional to the initial spatial deformation £j^, with the propor- 
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FIGURE 3. (color online) Scaled elliptic flow V2/ex as a function of (1 /S)dNch/dy (i.e. the charged 
multiplicity density per unit transverse overlap area S), for Pb+Pb and Au+Au collisions of different 
centralities at various collision energies, as compiled by the NA49 Collaboration |47]. The lines indicate 
the hydrodynamically predicted values at RHIC, SPS and AGS energies |4j|] (see text for details). 

tionality constant being controlled by the equation of state through the speed of sound 
c^. = dp/de [\]. The reason for this scaling is that the hydrodynamic equations of mo- 
tion are scale-free, and that therefore the initial condition for sets the scale for the 
final momentum anisotropy. This scaling assumes that the momentum anisotropy can 
fully develop, i.e. that the hydrodynamic evolution is not cut short by decoupling be- 
fore the spatial anisotropy Ex has completely disappeared. As discussed in the previous 
chapter, this assumption is justified near midrapidity in central and near-central Au-i-Au 
collisions at RHIC, but breaks down more and more severely as one goes to more periph- 
eral collisions, away from midrapidity, or to lower collision energies. This is at least part 
of the explanation of the difference between the data points in Fig. |3] and the horizontal 
lines labelled "HYDRO limits" [14]. 

These horizontal lines represent the fully developed hydrodynamic elliptic flow 
(scaled by the initial deformation e^) for minimum bias Au-i-Au or Pb+Pb collisions 
for three different effective sound speeds, reflecting the effective stiffness of our hadron 
gas -I- quark-gluon plasma EOS for collisions at (from right to left or bottom to top) 
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RHIC, SPS and AGS energies. They were extracted from the elliptic flow excitation 
function for b = 7fm Pb+Pb collisions given in Fig. 14 of [4], noting that b = lfm is 
roughly equal to the average impact parameter in minimum bias Au+Au collisions. The 
RHIC value for V2/ £x is lower than the ones for SPS and AGS energies since RHIC col- 
lisions spend more of the relevant time when the elliptic flow builds up in the soft phase 
transition region, compared to SPS and AGS where most or all of that time is spent in 
the stiffer hadron gas phase |4]. The increasing discrepancy between these lines and the 
data as one decreases the collision energy indicates the inability of the hadron gas phase 
(which gains more and more dynamical weight) to respond to the spatial deformation of 
the reaction zone. This is a late-stage effect and does not necessarily indicate a lack of 
early thermalization. 

However, some lack of early thermalization may still contribute to this discrepancy. 
Let us, for the sake of the argument, concentrate on Au-i-Au collisions at y/s = 200A GeV 
(solid red stars in Figure 01) and pretend that, in fact, the entire discrepancy between data 
and theory is due to inefficient thermalization at the beginning of the collision. Let us 
model this inefficiency in a somewhat black-and-white fashion, by assuming that the 
system evolves by free- streaming (i.e. without any interactions at all) for a time At, 
after which it instantaneously thermalizes and continues to evolve according to the laws 
of ideal fluid dynamics. As shown in [4] (for a slightly more general argument see EH]), 
free- streaming results in a radial growth of the source which reduces its initial spatial 
deformation according to 

e;,(To+AT) 

Due to the hydrodynamic scaling of the final V2 with the initial spatial deformation Sx, 
the final elliptic flow will be reduced by at least the same factor (and by even more 
if, after this delay, the hydrodynamic stage is no longer sufficiently long-lived to fully 
develop the elliptic flow). By assigning all of the "missing" elliptic flow in Figure |5] to 
this "delayed thermalization" effect, we can therefore extract a conservative upper limit 
for the thermalization time scale At. 

Note that even if (contrary to the above arguments based on Teaney's work lfl4ll ) our 
assumption were true that all of the "missing" elliptic flow in Figure |3] were due to 
inefficient thermalization in the initial partonic stage (rather than in the late hadronic 
stage), this procedure will still only yield an upper limit for the thermalization time At: 
In real life, thermalization doesn't happen suddenly after a period of no scattering at 
all, but gradually. Any type of scattering, however, will generate anisotropic evolution 
(albeit not as strongly anisotropic as in ideal fluid dynamics), i.e. faster expansion into 
the reaction plane than perpendicular to it, thereby reducing the spatial deformation 
(and thus the final V2) faster than given by Eq. ([T])- 

Let us now use this to extract an upper limit for the thermalization time At in ~ 7 fm 
Au-i-Au collisions (4th solid red star from the right in Fig.|3l). According to the Figure, 
for this impact parameter V2/£x is about 75% of the hydrodynamic value. Equating this 

1 /2 

with ex(To+AT)/ex(To), and inserting (r^)^,) =3.5 fm for b ^ 7 fm Au-i-Au, we find 
At ~ 2 fm/c for this upper limit. For second solid red star from the right, the experimental 
value corresponds to 90% of the hydrodynamic limit, and the resulting upper limit for 



1 + 



(cAt) 



2 



(1) 



the thermalization time is At ~ 1 fm/c. The right-most point in Figure El allows for a non- 
zero thermalization time only by virtue of its large error bar (and, of course, systematic 
uncertainties in the exact vertical position of the theoretical horizontal line). 

We see that, even if we assign all of the "missing" elliptic flow to inefficient early 
thermalization, the resulting window for the thermalization time scale is narrow and not 
larger than about 2 fm/c in minimum bias Au-i-A collisions and certainly less than 1 fm/c 
in central Au+Au collisions. This window becomes even narrower when accounting 
for the well-documented iffl] inefficiency of the late hadronic stage to generate elliptic 
flow. When taken together, these arguments yield very tight limits on the thermalization 
time scale in the dense matter formed early in the collision which, I think, are fairly 
represented (and not underestimated) by the value Teq = 0.6 fm/c used as starting time 
in our hydrodynamic simulation. (Let me remind you once again that this time was 
determined in a different way, using only information from central collisions, and that 
it therefore enters my network of arguments for early thermalization as an independent 
"data point".) 

It is, of course, unsatisfactory to see in Figure |3l the ideal hydrodynamic limit just 
being reached at one point, namely in the most central Au-i-Au collisions at the highest 
available collision energy. The natural reaction of any serious physicist should be that 
this single point coincidence must be accidental, and that the Figure really shows incom- 
patible tendencies in theory and data. The only problem with this "conservative" reaction 
is that it should be impossible to create more elliptic flow than predicted by ideal fluid 
dynamics! If at larger values for (1/5) dNch/dy the data for V2/£x indeed continued to 
grow significantly above the hydrodynamic limit (as suggested by the tendency of the 
data in Figure^), this would create very serious problems for theory. It is therefore very 
important to try to extend the range of ( 1 /S) dNch/ dy to larger values, by studying either 
semi-central Pb-i-Pb collisions at the LHC or central U-i-U collisions at RHIC [4, 46]. 
In order to support our present understanding of heavy-ion collision dynamics it will be 
essential to confirm that for ( 1/ S) dNch/dy > 30 the ratio vj/Sx saturates at the hydrody- 
namic limit, showing only a very weak growth between RHIC and LHC (resulting from 
a stiffening of the effective equation of state) |3] . 



Hadronization via quark coalescence - evidence for deconfinement 

We already mentioned that the right panel of Figure El indicates a characteristic 
difference between mesons and baryons in the way their elliptic flow breaks away at 
intermediate p± from the hydrodynamic behaviour at low p±. In the same intermediate 
p± region a similar meson-baryon split is seen in the nuclear enhancement factors, 
specifically the ratio Rcp between the various hadron yields (scaled by the number of 
binary nucleon-nucleon collisions) measured in central and peripheral Au-i-Au collisions 
isS]. It is quite striking that in this ratio the ^ meson data fall together with the much 
lighter kaons |53] and not with the almost equally heavy protons (s^ which instead 



follow the other measured baryons A, S, and Q. [ 53>.54] . This systematics indicates that 
hadron production at intermediate transverse momenta, 1 . 5 GeV /c<p±<5 — 6 GeV/c, 
is controlled by the number of valence quarks inside the hadron rather than by its mass. 



Mass scaling, as seen for the elliptic flow at low p±<l.5 GeV/c (see Fig.|21 right panel), 
is characteristic of thermalization and hydrodynamic flow; a scaling with the number 
of valence quarks, on the other hand, is not consistent with hydrodynamics and rather 
suggests a quark coalescense picture In fact, in the transition region around 1.5 — 
2.5 GeV/c where hydrodynamics begins to break down, a quark coalescence approach 
works extremely well, as I will now show. 

The coalescence model was previously developed for and applied to the formation 
of deuterons and other light nuclei from nucleons in low- and high-energy nuclear 
collisions (see fs^,"??] and references therein). The formalism developed there is most 
useful if the wave function of the cluster formed by the coalescence process has a narrow 
internal momentum distribution (i.e. it is relatively large). Clearly, hadrons have a much 
wider internal momentum distribution than deuterons, but if the total hadron momentum 
is sufficiently large (here we are interested in momenta > 1.5GeV/c), neglecting the 
internal relative momenta of the quarks inside the hadron may still be a reasonable 
approximation. This is how the coalescence model has bee n applie d to quark coalescence 



at RHIC by most authors ||58l|59|,|60] (for exceptions see 116 IL 16211 ). 

One starts from the following expressions for the invariant momentum spectrum for 
mesons M and baryons B [1561.15711 : 

dNM fp^da^{x) /■ 3 _M,a/3. . / , 1 \ r f 1 ^ 

^j3^ = Lj (2;rp /^^ V^P (g) fa{p + ^q,x) fp{p-^q,x), (2) 

dNs ^ fp^da^ix) f . . bm^y, ^ 

= ]^J^2nr-r '''''''' ^''^''^ 

X fa{p + \qi,x) fp {p + jq2.x) fY{p-j{qi+q2),x). (3) 



The i/a's are the internal wavefunctions of the hadrons, expressed as functions of the 
relative momenta of their valence quarks, and the /'s are the valence quark phase-space 
distributions. The sums extend over all spin, isospin and color channels contributing to 
the desired hadronic final state. The outer integral goes over the hadronization surface 
o{x). In the narrow wave function limit ^, = 0, and under certain assumptions about 
o{x) (e.g. that its spatial anisotropy is small so that the momentum anisotropy of the 
spectra is dominated by the momentum anisotropy of the quark distribution functions 
/(^ j8 jc)), these formulae imply that the azimuthal distributions of the hadrons are 
related to those of the quarks by 

dNM _ fdN^Y dNB _ [dNq~ 



d^p \d^pj d^p \d^p 
and that the elliptic flow coefficients (neglecting small nonlinear terms) satisfy ||6 

v^iP±) - m +vf (i^) , v'^iP^) ^ (^) +vP (i^) +vl (^) . (5) 
If all quark flavours carry the same elliptic flow, this simplifies further to ll60ll 

v'2(P±)^«v|(£^). (6) 



In other words, if we divide V2 for a hadron by the number n of its valence quarks and plot 
it against the scaled transverse momentum p±/n, we should obtain a universal function 
(^) describing the elliptic flow of the quarks just before hadronization. 
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FIGURE 4. (color online) Elliptic flow per valence quark, V2/n, as a function of transverse momentum 
per valence quark, p±/n, for pions, kaons, protons, A and S hyperons and their antiparticles. Data are 
from 200A GeV Au+Au collisions measured by PHENIX yjj (left panel) and STAR |49] (right panel). 
Universality of this curve suggests that it represents the partonic elliptic flow V2^'°"(p^'^°"), with no 
apparent difference between light and strange quark elliptic flow. 

Figure |4| shows that this universal scaling works excellently. There is no indication 
for different elliptic flows for light and strange quarks. The common parton elliptic 
flow shows the same generic form as the individual hadronic curves, following the 
hydrodynmic rise at low p± up to about p± = 150Me.Y where it breaks away from 
hydrodynamics, eventually saturating at a level of about 7% (for minimum bias Au+Au 
collisions) above = IGeV. Note that, on the parton level, hydrodynamics works only 
out to about twice the mean transverse momentum; the coalescence mechanism ensures 
that this is sufficient for mesons and baryons to behave hydrodynamically to much larger 
p^, namely about 1 .5 GeV for mesons and about 2.3 GeV for baryons. Pentaquarks ll63ll 
should therefore show hydrodynamic behaviour up to = 3.8 GeV/c, and their elliptic 
flow should saturate at 35%! 

Another important remark is that the extraction of the partonic flow from the rescaled 
hadronic elliptic flows is only possible because hydrodynamics breaks down at high p±. 
If the elliptic flow continued to follow the almost linear rise to higher p^, rescaling 
both vertical and horizontal axes would simply reproduce this linear curve, and nothing 
could be learned. Since breakdown of hydrodynamics produces a structure (shoulder) in 
the function V2{p±) and this shoulder moves under rescaling of the axes, the model can 
be tested and verified. The strongest evidence for the coalescence model comes from 
the fact that the shoulder of V2 happens for different collision centralities at different 
absolute values for V2 and at slightly different p± values; however, the quark number 
scaling of the shoulder position works at all impact parameters, i.e. different hadronic 
V2-curves collapse onto a single partonic V2 curve at each collision centrality ^4^. 

Careful readers will note that for <750MeV the extracted partonic V2 is not quite 



universal. In particular the differential partonic elliptic flow extracted from pions seems 
to be higher than that from the other hadrons (left panel in Figure |H). But this is only 
a reflection of the hydrodynamic mass splitting (remember: hydro works at low p±l) 
which, of course, cannot be scaled away by rescaling the axes. The coalescence model 
does not work at low p± (in particular, it violates entropy conservation because it 
reduces the number of entropy carriers by more than a factor 2!); it can only be used 
at sufficiently high p± where coalescence is a relatively rare process [60]. 

Clearly, the coalescence model makes essential use of the concept of deconfined 
(uncorrelated) quarks with independent phase-space distributions. I do not see how the 
fact that it works so perfectly can be interpreted in any other way than to provide at least 
indirect evidence for quark deconfinement in the early partonic phase. 



Conclusions: Have we seen the QGP? 

As the astute reader of these pages must have already concluded from my repeated 
use of the acronym QGP for the dense early collision stage, my own answer to this 
question is: YES! The quantitative success of the hydrodynamic model in describing 
the bulk of particles emitted from Au+Au collisions at RHIC provides compelling 
proof that thermalized matter at unprecedented energy densities (e) > 100 enm has been 
created, (^nm = 0. 125 GeV/fm^ is the energy density of cold nuclear matter inside atomic 
nuclei.) There is only one known viable theoretical concept that can describe thermalized 
matter at such energy densities: the Quark-Gluon Plasma. The mass splitting of the 
elliptic flow supports (but at the present accuracy does not prove) the hypothesis that 
this plasma is separated from normal hadronic matter by a phase transition (i.e. a soft 
region in the equation of state). The validity of the ideal fluid hydrodynamic model 
implies very severe constraints for the thermalization time scale and transport properties 
of this plasma: the required Teq < 1 fm/c cannot be obtained with standard perturbative 
QCD scattering mechanisms and requires (presently unknown) non-perturbative physics 
which makes the QGP created at RHIC behave like a strongly interacting, almost ideal 
fluid rather than a weakly interacting parton gas. 

One might worry that such non-perturbative processes might also modify the color 
deconfinement and chiral symmetry properties usually associated with the QGP, and 
transverse flow patterns do not, of course, provide any direct information about these 
characteristics. I showed that indirect arguments for color deconfinement can be ex- 
tracted from the pattern in which the hydrodynamic description of hadron spectra begins 
to break down at intermediate transverse momenta: This pattern is consistently described 
by a quark-coalescence model, in which collective flow patterns are first imprinted on the 
momentum distributions of deconfined (uncorrelated) quarks in the early partonic phase 
and then transferred to the late hadronic phase by a quark recombination mechanism. 

Another indirect argument going in the same direction is based on the success of the 
grand canonical thermodynamic ensemble in describing all the measured hadron abun- 
dance ratios L30.1 . This approach works extremely well even for multistrange hyperons 
and antihyperons, with a fully equilibrated strange particle sector (in contrast to ele- 
mentary particle collisions and nuclear collisions at lower energies where strangeness is 



found to be suppressed relative to the grand canonical approach and requires taking into 
account, via canonical constraints and/or additional strangeness suppression factors, lo- 
cal correlations between strange quarks and antiquarks |64]). This implies that, say, in 
the production of an Q. hyperon, strangeness can be balanced by kaons or anti-hyperons 
anywhere in the fireball volume, including its far opposite edge. In other words, any 
memory of the fact that the microscopic strangeness production process in QCD is lo- 
cal, i.e. s and s are created pairwise in a point, has been completely lost by the time 
strange hadrons form, and at hadronization the strange quarks and antiquarks are com- 
pletely uncorrected and equally likely to be found at any point in the fireball. It is very 
difficult, if not impossible, to explain this without invoking strong rescattering dynamics 
among deconfined quarks and gluons before the onset of hadron formation. 
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